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Abstract 

The nature of decoupling in the mixed phase of extremely type-II layered superconductors 
is studied theoretically through a duality transformation of the layered XY model with 
frustration. In the limit of weak coupling, we generally find that the Josephson effect is 
absent if and only if the phase correlations within isolated layers are short range. In the case 
specific to uniform frustration, we notably identify a decoupled pancake vortex liquid phase 
that is bounded by first-order and second-order decoupling lines in the magnetic field vs. 
temperature plane. These transitions potentially account for the flux-lattice melting and 
for the flux-lattice depinning that is observed in clean high-temperature superconductors. 
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Introduction. High-T c superconductors are perhaps the best known example of an ex- 
tremely type-II layered superconductor. [1] The vortex lattice that exists in the mixed phase 
of clean oxide superconductors melts through a first-order transition for magnetic fields 
applied perpendicular to the layers. [2] The first-order line, H± = H m (T), begins at the 
zero-field critical point, T c , but it ends strangely in the middle of the phase diagram. The 
depinning line T = Td p (H±), which marks the point at which the flux lattice depins itself 
through thermal excitations, appears to be independent of this melting line. [3,4] The first- 
order transition is commonly interpreted as either a vortex lattice melting transitional] 
or as a layer decoupling transition. [5,6] The respective theoretical approaches have been 
primarily based on elastic medium descriptions of the vortex matter that notably omit 
topological defect excitations. [7,8] These theories also fail to account for the multicritical 
point mentioned above in the absence of disorder. [9] 

In this paper, we shall model the thermodynamics of the interior of the mixed phase 
in extremely type-II anisotropic superconductors by the layered XY model with uniform 
frustration. [10-12] The nature of the Josephson effect between layers is uncovered through 
a partial duality analysis of the XY model. [13-17] A subsequent weak-coupling analysis at 
high perpendicular fields yields a low-temperature phase made up of coupled 2D vortex 
lattices that is separated from a decoupled liquid of planar vortices at high temperatures 
by a second order melting line. It is further argued that the decoupled liquid phase 
experiences a first-order phase transition into a coupled solid phase as the perpendicular 
field is lowered into the strong-coupling regime. [16,17] The phase diagram that results (see 
Fig. 1) is compared with recent experimental reports of multicritical phenomena in the 
mixed phase of high-temperature superconductors. [3,4] 

Layered XY Model. Before considering the uniformly frustrated case that describes 
the mixed phase of extremely type-II superconductors, let us first perform a duality analysis 
of the layered XY model without frustration. [18] The corresponding energy functional 
reads 
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where <f>(r, I) is the superconducting phase at a point r in layer /. Above , A M denotes the 
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nearest-neighbor difference operator along the ft direction, while represents a purely 
longitudinal vector potential. The Josephson penetration length j'a provides a natural 
scale for the model (1) in the limit of weak inter-layer coupling, in which case the model 
anisotropy parameter 7' = (J\\/ J±) 1 ^ 2 is much larger than unity. Here, a denotes the 
(square) lattice constant for each layer. Last, any generalized phase auto-correlation func- 
tion set by an integer source field, p(r), is related to the corresponding partition function 

Z f Y [p] = J V ( pe- E xl-/^ T e^ p<t> (2) 

by the quotient 
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We shall now employ the well-known dual representation[13] of the XY model (1) based 
on the Fourier series expansion 

e p C ose = ^^° = _ oo /| n |(/3)e ine of the Gibbs distribut ion in 
terms of modified Bessel functions, I n (x). This identity allows the phase variables to be 
integrated out of (2). Resummation over the parallel link fields n x and n y then yields the 
form 

z&lp] = ifm ■ z cg \p\ ■ n^z&io] (4) 

for the partition function of the layered XY model in terms of a product of a layered 
Coulomb gas ensemble (CGE) 

z cg [p]= £ v^-^tiCM-e-^^ (5 ) 
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with N isolated XY model layers (J± = 0).[17] Above, n z (r) is the remaining integer link 
field between adjacent layers that is effectively restricted to take values n z (r) = 0, ±1 in 
the weak-coupling limit Jj_ ksT. Indeed, each configuration in the CGE (5) is weighted 
by the bare fugacity yo = J±/2ksT raised to the power N[n z ] = J2ri \ n z(r,l)\ equal to 
the total number of n 2 -charges. The latter so-called fluxons physically represent vortex 
rings in between adjacent layers. [14] Each configuration is further weighted by a product 
of generalized auto-correlation functions Ci[q{\ = (exp i^^qi{r)(f){r,l) ) over each 

\ L J / j^=o 

layer / in isolation that are evaluated at sources qiif) = p(f, I) + n z (r, I — 1) — n z (r, /). Last, 
we have the parameters /3^±_ = J||,_l//cb^ 5 while M' denotes the total number of rungs 
between layers. 



To proceed further, it is instructive to consider a single neutral pair of unit n z charges 
that lie in between layers V and V + 1 in the absence of an external source, p = 0, with 
the negative and positive charges located at planar sites r*i and f 2 , respectively. The 
gauge-invariant product over intra-layer autocorrelation functions in the layered CGE (5) 
then reduces to the product \Ci'(fi — F2) | ■ |CV+i(fi — f 2 )\ of the corresponding phase 
autocorrelation functions, 

C/(rL,r 2 ) = ^exp z0(ri,Z) - i(p(f 2 J) ^ _^ , (6) 

within isolated layers. This function varies algebraicly with the separation as 

\Ci(f)\=go(r /\r\r 2D for |f| « £ vx , (7) 

while it decays exponentially with the separation as 

|Ci(r)l=0oexp(-|f|/£vx) for |r=1 » £ vx . (8) 

Here, t] 2 d = Vsw + Vv* is the 2D correlation exponent inside layer /, where rj sw = (2tt(3\\)~ 1 
and rj vx are the respective spin-wave and vortex contributions. Also, £ vx denotes the 2D 
phase correlation length, while the length r = a/(2 3 / 2 e 7 ) is set by Euler's constant, 7. 
The effective layered CGE (5) therefore takes the form 
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in the limit of dilute fluxon (n z ) charges. At separations within a correlation length, 
\ri — r 2 \ (vx, the fluxons experience a pure Coulomb interaction (V^} =0). At large 
separations \r\ — r 2 \ >> ^ vx , on the other hand, the fluxons experience a pure (j] 2 d = 0) 
confining interaction (f\, r 2 ) = \ f\ — r 2 1 / between those points fi and r 2 in layer 

/ that are connected by a string [see Eq. (8) and ref. [19]]. Below, we determine the 
thermodynamic nature of the superconducting and "normal" phases that correspond to 
Eqs. (7) and (8), respectively. 

Consider first the case where (quasi) long-range intra-layer phase correlations are 
present: £ vx = 00 and V^} = 0. Summing independently over charge configurations of 



the CGE (9) that are restricted to take values n z = 0, ±1 at each site plus an appropri- 
ate Hubbard-Stratonovich transformation[12] yields an equivalent renormalized Lawrence- 
Doniach (LD) model[l] set by the continuum energy functional 

r r N i Ar_1 
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Here, Jy = ksT /2txt]2d is the macroscopic phase rigidity of an isolated layer, [20] while 
Aq = a(J\\/goJ±) 1 ^ 2 is the renormalized Josephson scale. The above continuum description 
(10) is understood to have an ultra-violet cut-off on the order of the lattice constant, 
ro ~ a. Eq. (10) is known to sustain a macroscopic Josephson effect at temperatures below 
hsT* — 4-7T J||.[14,15] Also, the factorization (4) plus the continuum limit (10) indicate the 
expression 

(cos0 w )-^ + / o ^j (11) 

for the local Josephson coupling, where (r) = <j>(f, l + l) — 0(r, /) — A z (f) is the gauge- 
invariant phase difference between consecutive layers. An analysis of the double-layer case 
yields an effective anisotropy scale Aj = Aq/2 1 / 2 and the limiting values r\ — > r)2D and 
fo — * 9o for the exponent and for the prefactor at low temperatures, i]2D *C 1 (see ref. 
[17]). Expression (11) indicates that a crossover to strong coupling, (cosc^+i) ~ 1, must 
therefore take place in the isotropic regime a ~ A and/or at a temperature of order the 
Josephson scale Tj = J±/ks- 

Consider next the case in which intra-layer correlations are short range: £ vx < oo. 
Inter-layer fluxon (n z ) pairs are then bound by a confining string (8). Application of the 
CGE (5) yields Koshelev's formula[ll] 

(e^) =y J d 2 r\d(r)\ ■ \C l+1 (r)\/a 2 (12) 

for the local Josephson coupling in this decoupled phase (see ref. [17]). Substitution of the 
form (8) for the short-range intra-layer auto correlator in turn yields the explicit formula 

(cos^ +1 ) = ^^ P ± (13) 

for the local Josephson coupling. It can be shown that the macroscopic Josephson effect 
is absent in this decoupled phase. [16,17] 



We can now determine the phase diagram of the layered XY model. [18] It is well known 
that an isolated XY layer looses phase coherence above a temperature /c^Ti 21 ^ = § «/||-[13] 
The transition is driven by the unbinding of vortex/anti- vortex pairs and it is second- 
order. [21] Consider first the weak-coupling limit (cosc^z+i) — > 0. By the previous analysis, 
we then conclude that the layers are Josephson coupled at low temperature T < Tc 2D ^ 
following the renormalized LD model (10), while they are decoupled at high temperature 
T > Tc 2D \ Eq. (13) indicates, also, that the selective high-temperature expansion breaks 
down ((cos (pi,i+i) ~ 1) in the decoupled phase at a temperature T x set by the identification 
of length scales A ~ £ vx (r x ). The second-order transition therefore takes place at a 
temperature T c that lies inside of the dimensional crossover window T C (2D) < T < T x 
for large yet finite anisotropies 7' ^> 1. These conclusions agree with what is presently 
understood for the layered XY model. [18] 

Uniformly Frustrated Case. Consider now the layered XY model (1) in the presence 
of a uniform frustration, = (0, b±x, 0), which describes the mixed phase of an extremely 
type-II superconductor in a field B± = (<& /2na)b± aligned perpendicular to the layers. [10- 
12] The duality analysis just performed can then be repeated wholesale, yet with the 
following modifications: [16,17] (a) Each isolated XY layer now undergoes a 2D melting 
transition that is mediated by the unbinding of dislocation pairs[21] at a temperature [22] 
ksTm^ — ^||/20, while (b) the ultra-violet cut-off, tq, of the renormalized LD model (10) 
is now of order the average spacing between planar vortices, a vx = (&o/B±) 1 / 2 . (The CGE 
[Eqs. (7)-(9)] must be coarse-grained up to the new ultra-violet scale, a vx .) In the weak- 
coupling limit (cos0^ + i) — > reached at high fields B± 3> B\_ = $ Mo> we therefore 
have iV 2D vortex lattices that show a macroscopic Josephson effect at low temperatures 
T < Tm , while a decoupled liquid of planar vortices exists at high temperatures T > 
Tm D \ At large yet finite anisotropies, 7' 3> 1, it is useful once again to determine the 
temperature scale T x at which point the selective high-temperature expansion (13) for the 
local Josephson coupling breaks down. This again takes place roughly when A and £ vx 
are comparable. As in the previous case without frustration, we therefore expect a second- 
order melting transition at a temperature T m that lies inside the dimensional crossover 
window[18] T% D) <T <T x . 



Yet what happens when the local Josephson coupling (cos <t>i,i+i) approaches unity as 
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the field B± is lowered? It is useful to first define a decoupling contour 

(cos0 M+ i) = (cos0 M+ i)£>, (14) 

in the T-B± plane. Numerical simulations indicate that (cosc^z+i) d is a constant less than 
but of order unity. [11] The result (11) for the local Josephson coupling in the coupled phase 
yields (i) a contour line at temperatures of order the Josephson energy, ksTj = Jj_, for high 
perpendicular fields B± ^> B^, and (ii) a contour line at perpendicular fields of order B*^ 
for temperatures near Tm D K Since no phase transition is possible in the screened CGE (9) 
that describes the coupled phase, this contour line must therefore represent a crossover into 
a flux-line lattice regime that exists at lower temperatures and fields. At high temperatures 
T > T x inside the weak-coupling regime of the decoupled phase, the string interaction (8) 
binds together fluxon-antifluxon pairs into stable dipoles of dimension £ vx that do not 
overlap. It can be shown that such fluxon pairs begin to overlap and dissociate in the 
vicinity of the decoupling contour (14) for temperatures that lie outside of the 2D critical 
regime (£ vx ~ a vx ).[17] We therefore expect a first-order transition along this line due to 
the absence of a nearby divergent length scale. 

The above results are summarized by the schematic phase diagram shown in Fig. 1. 
The phenomenology J± oc (T c — T)/T c for the Josephson energy in the vicinity of the zero- 
field transition at T c yields the same linear temperature dependence for the first-order 
decoupling field, [2,5,6] Ho(T) [see Eqs. (13) and (14)]. Eq. (11) also implies that the 
decoupling contour (14) continues into the coupled phase at temperatures T < T m , there 
representing a crossover. The first-order line (14) must therefore end at the second-order 
melting line, T = Tm. The melting curve then continues up in field along the latter second- 
order line. Finally, it is possible that a vestige of the second-order melting transition at 
high fields B± > B\ persists down into the low-field region in the form of a crossover (see 
Fig. 1). 

Discussion. The first-order flux-lattice melting line in clean high-temperature su- 
perconductors is experimentally observed to end at a multi-critical point. [2,3] The phase 
diagram proposed in Fig. 1 for the mixed phase of extremely type-II layered superconduc- 
tors also exhibits a multicritical point at a temperature and field T ~ T^^ and B ~ B*^. 
This point agrees qualitatively with the high-T c phase diagram. [2,5] Also, bulk pinning is 
observed to become relatively strong in clean oxide superconductors at low temperatures 
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and high fields, T < T£ u> and H± > B* ± .[3] The quasi-2D vortex lattice phase that is 
identified in Fig. 1 coincides with this regime. It is a smectic (super) solid, [8] and it can 
thus adjust better to a random landscape of point pins than a rigid vortex lattice. Recent 
observations of muon spin resonance in the mixed phase of clean high-temperature super- 
conductors are consistent with such a picture of enhanced point pinning due to dimensional 
crossover. [23] We remind the reader that all effects due to magnetic screening[24] have been 
completely neglected in the present theory. 

In conclusion, the phase diagram proposed in Fig. 1 for the layered XY model with 
uniform frustration is strikingly similar to that of the mixed phase in high-T c supercon- 
ductors. [2,3] It is perhaps more important to point out, however, that the above duality 
analysis yields only two thermodynamic phases at weak coupling: a coupled superconduc- 
tor (7) and a decoupled "normal" state (8). This indicates that neither the Friedel scenario 
(decoupled superconducting layers) nor the flux-line liquid state (coupled normal layers) 
are thermodynamically possible in the absence of disorder. [1] 

The author thanks M. Dodgson for correspondence and B. Rosenstein for discussions. 



8 



References 



1. G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. Larkin, and V.M. Vinokur, Rev. 
Mod. Phys. 66, 1125 (1994). 

2. E. Zeldov, D. Majer, M. Konczykowski, V.B. Geshkenbein, V.M. Vinokur, and H. 
Shtrikman, Nature 375, 373 (1995). 

3. D.T. Fuchs, E. Zeldov, T. Tamegai, S. Ooi, M. Rappaport and H. Shtrikman, Phys. 
Rev. Lett. 80, 4971 (1998); B. Khaykovich, E. Zeldov, D. Majer, T.W. Li, PH. Kes, 
and M. Konczykowski, Phys. Rev. Lett. 76, 2555 (1996). 

4. T. Shibauchi, T. Nakano, M. Sato, T. Kisu, N. Kameda, N. Okuda, S. Ooi, and T. 
Tamegai, Phys. Rev. Lett. 83, 1010 (1999). 

5. L.I. Glazman and A.E. Koshelev, Phys. Rev. B 43, 2835 (1991). 

6. L.L. Daemen, L.N. Bulaevskii, M.P Maley and J.Y. Coulter, Phys. Rev. B 47, 11291 
(1993). 

7. M. Feigel'man, V.B. Geshkenbein, and A.I. Larkin, Physica C 167, 177 (1990). 

8. E. Frey, D.R. Nelson, and D.S. Fisher, Phys. Rev. B 49, 9723 (1994). 

9. B. Horovitz and T.R. Goldin, Phys. Rev. Lett. 80, 1734 (1998). 

10. Y.H. Li and S. Teitel, Phys. Rev. B 47, 359 (1993); ibid 49, 4136 (1994). 

11. A.E. Koshelev, Phys. Rev. Lett. 77, 3901 (1996); Phys. Rev. B 56, 11201 (1997). 

12. J. P. Rodriguez, J. Phys. Cond. Matter 9, 5117 (1997). 

13. J.V. Jose, L.P Kadanoff, S. Kirkpatrick and D.R. Nelson, Phys. Rev. B 16, 1217 
(1977). 

14. S.E. Korshunov, Europhys. Lett. 11, 757 (1990). 

15. J. P. Rodriguez, Europhys. Lett. 31, 479 (1995). 

16. J. P. Rodriguez, "Nature of Decoupling in the Mixed Phase of Extremely Type-II 
Layered Superconductors" (|cond-mat/9906199| ) . 

17. J. P. Rodriguez, "Theory for Decoupling in the Mixed Phase of Extremely Type-II 
Layered Superconductors", ( |cond-mat / 0006 176] ) . 

18. S. Hikami and T. Tsuneto, Prog. Theor. Phys. 63, 387 (1980); see also S.R. Shenoy 
and B. Chattopadhyay, Phys. Rev. B 51, 9129 (1995). 

19. A. Polyakov, Phys. Lett. 72 B, 477 (1978). 

20. D.R. Nelson and J.M. Kosterlitz, Phys. Rev. Lett. 39, 1201 (1977). 

9 



21. J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 (1973). 

22. S.A. Hattel and J.M. Wheatley, Phys. Rev. B 51, 11951 (1995). 

23. J.E. Sonier, J.H. Brewer, R.F. Kiefl, D.A. Bonn, J. Chakhalian, S.R. Dunsiger, W.N. 
Hardy, R. Liang, W.A. MacFarlane, R.I. Miller, D.R. Noakes, T.M. Riseman and C.E. 
Stronach ( |cond-mat/9906449|) . 

24. G. Blatter, V. Geshkenbein, A. Larkin and H. Nordborg, Phys. Rev. B 54, 72 (1996). 



10 



Figure Caption 



1. Shown is a schematic phase diagram for the uniformly frustrated XY model (1) made 
up of a finite number of weakly-coupled layers. The cross-over transition that is 
represented by the horizontal dashed line at low temperatures Tj < T < is 
broad (see ref. [17]). 
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